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ANALYSIS AND FRACTIONATION OF PARTICLES NEAR SURFACES 
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5 Application Serial No. 09/813,571 filed March 21, 2001 , which is incorporated 
h e rein by r e f e renc e in its e ntirety . 
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FIELD OF THE INVENTION 

The present invention generally relates to the fields of analytical chemistry 
and materials characterization. Specially, the present invention relates to 
1 5 fractionation of heterogeneous mixtures of particles on the basis of characteristic 
frequency dependence of particle polarization and also to analysis of 
electrochemical properties of particles. 

BACKGROUND OF THE INVENTION 

LI Methods of Fractionation and Sorting 

20 Among analytical separation techniques, methods of flow cytometry have 

been used in sorting of particles (including cells) within a given heterogeneous 
population. These methods, however, require a dedicated instrument capable of 
selecting and physically separating from the population those particles which satisfy 
a given selection criterion. For example, the objective of cell sorting by flow 

25 cytometry usually is the selection of those cells within a population displaying a 
characteristic fluorescence signal. 

1.2 Methods of Particle Analysis 

Particle analysis represents a standard procedure of the analytical chemical 
repertoire that is used to determine physico-chemical and compositional properties 
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[Hunter, "Introduction to Modern Colloid Science", Oxford University Press, 
Oxford, UK, 1993]. The extensive repertoire of techniques reflects the ubiquitous 
use of particles of a wide range of sizes, shapes, composition and chemical 
reactivity in many scientific and industrial applications ranging from medical 
5 diagnostics to cosmetics. Characterization of particles is useful in guiding and 
optimizing production as well as chemical modification particularly of surface 
properties that determine the stability of particle suspensions and the interaction of 
particles with molecules in the surrounding medium. 

Techniques of the standard repertoire may be grouped as follows: first, 

10 sedimentation and centrifiigation methods, electroacoustics, light scattering, 
hydrodynamic methods and dielectric spectroscopy; second, electrical pulse 
counting, flow cytometry and electrophoretic zeta potential measurements; third, 
specialized methods including dielectrophoresis (DEP), including its combination 
with field-flow fractionation methodology. 

1 5 Methods in the first group apply to bulk suspensions and therefore deliver an 

average of the measured quantities for a large number of particles. While desirable 
in situations in which large industrial batches are to be characterized, bulk 
measurements generally do not lend themselves to miniaturization and are not well 
suited to the characterization of small numbers of particles, particularly when 

20 particle-to-particle variations are of interest or when multiple particle parameters are 
to be determined simultaneously. 

Methods in the second group use a "single- file" serial format of 
measurement. While these methods determine the properties of individual particles, 
their implementation, in order to attain sufficiently high processing rates, requires 

25 high flow rates of a carrier gas or fluid and high-speed data recording and read-out 
electronics, generally rendering the equipment complex and expensive. This is so 
especially when multiple parameters are to be determined for each particle. For 
example, multi-color detection in flow cytometers can require the use of multiple 
lasers and multiple photomultiplier tubes. The determination of the zeta potential 

30 along with particle size in state-of-the-art equipment requires the measurement of 
electrophoretic mobility of particles traveling along a narrow capillary in response 
to a DC voltage applied along the channel as well as the application of light 
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scattering with associated light source and read-out, as in the ZetaSizer (Malvern 
Instruments, Southborough, MA). Miniaturization, as in the case of flow cytometry 
[Fu, AX, Spence, C, Scherer, A., Arnold, F.H. and Quake, S.R. A Microfabricated 
Fluorescent- Activated Cell Sorter, Nature Biotechnology, Vol. 17, November 1999, 
5 1 109-1 1 1 1], retains the serial mode of processing and thus encounters the constraint 
of limited throughput. 

Methods in the third group include various dielectrophoretic techniques to 
characterize, classify and fractionate low conductivity particle suspensions [Becker, 
F.F., Gascoyne, P.R.C., Huang, Y. and Wang, X-B, Method and Apparatus for 

10 Fractionation Using Generalized Dielectrophoresis and Field Flow Fractionation, 
US Patent No. 5,888,370, March 30, 1999; Rousselet, J., Markx, G.H. and Pethig, 
R. , Separation of Erythrocytes and Latex Beads by Dielectrophoretic Levitation and 
Hyperlayer Field-Flow Fractionation, Colloids and Surfaces A: Physicochemical 
and Engineering Aspects, 140 (1998) 209-216; Pethig, R., Markx, G.H., Apparatus 

1 5 for Separating by Dielectrophoresis, US Patent No. 5,8 1 4,200, Sep 29 th , 1 998; Pohl, 
H.A., Continuous Dielectrophoretic Cell Classification Method, US Patent No. 
4,326,934, April 27 th 1982; Parton, A., Huang, Y., Wang, X-B., Pethig, R., 
MacGregor, A.R., and Pollard-Knight, D.V., Methods of Analysis/Separation, US 
Patent No. 5,653,859, Aug. 5, 1997; Crane, S., Dielectrophoretic Cell Stream Sorter, 

20 US Patent No. 5,489,506, Feb. 6 th , 1996; Benecke, W., Wagner, B., Hagedorn, R., 
Fuhr, G., Muller, T., Method of Continuously Separating Mixtures of Microscopic 
Dielectric Particles and Apparatus for carrying through this method, US Patent No. 
5,454,472, Oct. 3 th 1995]. All publications which are cited throughout the 
application are hereby incorporated by reference in their entirety. 

25 Coupled with a suitable imaging system, techniques within this group can 

provide single particle information. The particle classification is achieved by 
exposing the particle mixture to a non-uniform AC electric field which is generated 
by applying an AC voltage to multiple sets of planar patterned microelectrodes. 
Such a scheme in general requires a complex electrode signaling/addressing 

30 scheme. Also, a continuous mode of operation requires, in the simplest case, 
coupling to an external flow or the use of field flow fractionation. Simple 
dielectrophoretic setups do not generally allow for quantitative evaluation for the 
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surface potential or surface conductivity of particles. While these quantities maybe 
determined from electrorotation spectra of single particles, this approach requires a 
complex and time-consuming set-up to confine a single particle at the origin of a 
rotating electric field whose frequency is scanned while recording the rotational 
5 motion of the confined particle. Aside from the complexity of the equipment 

required to implement the requisite experimental configuration, the low throughput 
of the method presents a serious disadvantage when averages of multiple particles 
are desired. 

SUMMARY OF THE INVENTION 

1 0 The present invention provides a method and apparatus for the fractionation 

of a mixture of particles and for the analysis of particles, based in part on the 
technology designated "LEAPS", described inPCT International Application 
No.WO 97/40385 to Seul. The present invention uses LEAPS (referring to "Light- 
Controlled Electrokinetic Assembly of Particles near Surfaces") to fractionate a 

1 5 heterogeneous mixture of particles on the basis of, e.g.: morphological properties 
including size and shape; surface composition relating to certain electrochemical 
properties, including surface potential and surface conductivity, and specific and 
non-specific adhesion to exposed surfaces provided within the apparatus. As used 
herein, the term "fractionation" refers to separation or sorting of particles, including 

20 separation of multi-component particle mixtures into their constituent 

subpopulations. Particles may be separated on the basis of the characteristic 
frequency dependence of the particle polarization, said dependence reflecting 
physical or chemical properties of the particles. 

The present invention also provides for the use of LEAPS to determine 

25 various physical or chemical properties of particles, including morphological 

properties, surface electrochemical properties and surface chemical properties such 
as the particles adhesiveness for given surfaces of the apparatus or for other 
particles. 

The present invention provides a method and apparatus for fractionation of 
30 particle mixtures on the basis of the respective characteristic relaxation frequencies 
of the constituent particles or particle populations. In certain embodiments of such 
method, a plurality of particles are suspended at an interface between an electrolyte 
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solution and a light-sensitive electrode. The particles comprise at least two types, 
each type having a distinguishable relaxation frequency. An electric field, having a 
frequency that is less than or equal to the relaxation frequency of at least one of said 
particle types but greater than the relaxation frequencies of other particle types, is 
5 generated at the interface by application of an AC voltage. In addition, the interface 
is illuminated with a predetermined light pattern. The illumination in combination 
with the electric field produces fractionation of the particles having relaxation 
frequencies greater than or equal to the frequency of the electric field from other 
particles, e.g., by assembly of particles having a relaxation frequency greater than 

10 the frequency of the applied electric field resulting in formation of a planar array of 
substantially one layer of particles in an area on the electrode designated by the 
pattern of illumination. In certain embodiments, the illumination affects the 
electrolyte-electrode interfacial impedance. 

In certain preferred embodiments, the electrode comprises a silicon electrode 

1 5 which is coated with a dielectric layer. In certain preferred embodiments, an 
additional electrode is provided such that the light-sensitive electrode and the 
additional electrode are substantially planar and parallel to one another and 
separated by a gap (e.g., in a sandwich configuration), with the electrolyte solution 
containing the particles being located in the gap. 

20 In certain other embodiments of the present invention, patterned electrode is 

used instead of or in combination with the illumination to fractionate a mixture of 
particles. In one such example, a first electrode positioned in the first plane and a 
second electrode positioned in the second plane different from the first plane is 
provided with a gap between the electrodes, an arrangement herein referred to as a 

25 sandwich configuration. A mixture of particles is suspended in an electrolyte 

solution and located in the gap. The second electrode comprises a planar electrode 
having a surface and an interior, either of which may be patterned to modify the 
spatial distribution of the interfacial (interface between the electrolyte solution and 
the second electrode) electric field. When an electric field is generated between the 

30 first and the second electrode by applying an AC voltage between the two, the 
electric field in combination with the patterning of the electrode produces 
fractionation of the particles according to relaxation frequency; for example, 
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particles having a relaxation frequency greater than the frequency of the applied 
electric field move toward designated regions of the second electrode surface as 
provided by illumination or patterning and therein form a planar array of 
substantially one layer. In preferred embodiments, the second electrode may be 
5 patterned by spatially modulated oxide growth, surface chemical patterning or 
surface profiling (which is described in U.S. Serial No. 09/171,550, filed October 
26, 1998, now U.S. Patent 6,251,691 issued June 26, 2001); preferably, the 
patterning affects the impedance of the electrolyte solution and second electrode 
interface. In preferred embodiments, the second electrode comprises a silicon 

10 electrode which is coated with a dielectric layer. 

In certain embodiments, the difference in relaxation frequency of particles 
within a heterogeneous mixture corresponds to a difference in particle size or 
surface composition, thus permitting fractionation of the mixture into homogeneous 
subpopulations of particles according to size or surface composition, respectively. 

1 5 The process of fractionation may also be monitored using a video detector or 
camera. 

Accordingly, the present invention provides a method and apparatus for 
fractionation of particle mixtures based on the characteristic relaxation frequencies 
of said mixture's constituent particles or particle subpopulations. The present 

20 invention provides for the determination of the frequency dependence of the 
electrohydrodynamic forces that govern particle transport into designated low- 
impedance areas of the electrode and mediate formation of particle arrays in those 
designated areas. As described herein, a particle relaxation frequency, co R , is a 
property that is associated with the particle-solution interface and may be 

25 determined from the condition that array assembly requires that the frequency 
©<o)r,: particles not meeting this criterion will be excluded from low-impedance 
areas. 

The present invention also provides for a method and apparatus for the 
analysis of particles in order to determine various physical or chemical properties of 
30 said particles, including size, morphology, surface potential and surface 

conductivity. For example, the relaxation frequency and/or maximal velocity of 
transport of the particles are used for particle characterization. 
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In certain embodiments of the present invention, a plurality of particles 
suspended in an interface between an electrolyte solution and a light-sensitive or a 
patterned electrode are provided, as described above in connection with particle 
fractionation. When an electric field is generated at the interface by the application 
5 of an AC voltage, the frequency of the electric field is adjusted to produce formation 
of a planar array of substantially one layer of particles in an area designated by the 
pattern of illumination or electrode patterning. From the analysis of particle 
transport in the course of array assembly, the relaxation frequency and the maximal 
velocity of transport of the particles may be determined. Alternatively, the particle 

1 0 relaxation frequency also may be determined by analyzing the response of the array 
configuration to changes in the frequency of the applied electric field. 

The present invention also provides for the determination of the frequency- 
dependent and voltage-dependent maximal velocity, Vmax attained by particles 
crossing impedance gradients in the course of array assembly. As described herein, 

1 5 Vmax may be determined by image analysis and particle tracking. In certain 
embodiments, the combination of cor and V max provides for the simultaneous 
determination of the particles' surface ("zeta") potential, <|> s , and surface 
conductivity, a, e.g., based on numerical analysis of phenomeno logical equations of 
particle motion described herein. 

20 Once these values are found, they may be used to characterize the particles, 

e.g., determine the zeta potential of the particles and/or the mobility of ions or 
molecules within the electrolyte solution in proximity to the particle surface. In 
preferred embodiments, the relaxation frequency and the maximal velocity are used 
together in particle characterization. In certain preferred embodiments, the present 

25 invention permits the simultaneous determination of multiple particle properties, 
such as zeta potential and ionic or molecular mobility. In other embodiments, the 
present invention permits the determination of one particle property of interest 
independently from the determination of others. For example, in certain 
embodiments, the ionic or molecular mobility associated with the particle-solution 

30 interface may be determined from the relaxation frequency alone. 

The term "particles," as used herein include colloidal beads, eukaryotic and 
procaryotic cells, micelles, vesicles, and emulsion droplets. In preferred 
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embodiments, the particles comprise colloidal beads, or eukaryotic or procaryotic 
cells. 

In contrast to prior art methods of fractionation of particle mixtures such as 
field flow fractionation, the methods of the present invention provide a parallel 
5 format based on the manipulation of a multiplicity of individually imaged particles 
in a simple sandwich geometry which provides for longitudinal and transverse 
electro hydrodynamic forces to mediate fractionation on the basis of the frequency- 
dependent polarization of the particles of interest. Specifically, the methods of the 
present invention provide for the fractionation of one subpopulation of particles at a 

1 0 time from the remainder. 

In contrast to prior art methods of sorting such as fluorescence-activated 
particle and cell sorting, the methods of the present invention provide for a parallel 
process accommodating a multiplicity of particles thereby enhancing processing 
speed. The methods of present invention also are well suited for miniaturization. 

15 In contrast to prior art methods of particle analysis which are applied to bulk 

particle suspensions and yield ensemble averages of the particle properties of 
interest, the methods and apparatus of the present invention provide for the 
miniaturization of the processes of analysis. In contrast to prior art serial methods of 
analysis of particle mobility, the methods of present invention invoke a collective 

20 phenomenon of particle array assembly as the basis for a parallel mode of 

simultaneous analysis of multiple particle properties accessible either to imaging 
and image analysis or to light scattering. The methods of the present invention 
provide for the simultaneous, yet independent determination of morphological 
parameters including size, orientation and polydispersity, as well as for 

25 electrochemical parameters including surface potential and surface conductance. 

In contrast to these prior art methods, the implementation of the fractionation 
methods of the present invention are readily suitable for miniaturization and 
integration with additional preparative and analytical procedures in an integrated 
analytical environment. 

30 BRIEF DESCRIPTION OF DRAWINGS 

Fig. 1 is an illustration of the response of a polarizable double layer around a 
spherical particle to a time-varying electric field. The progression of polarization 
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states occurring in a single cycle is shown in panels (a) - (e). The shading indicates 
the distribution of the predominant type of ion in the diffuse layer. 

Fig. 2 is an illustration of the variation of the polarization factor of a particle 
with frequency, as evaluated numerically on the basis of the theoretical model 
5 described herein and assuming the indicated values of the particle surface potential, 
a particle diameter of 2.8 \im and a solution of ionic strength 0.1 mM. 

Fig. 3(A) is an illustration of the dependence of the relaxation frequency on 
surface potential, as evaluated numerically on the basis of the theoretical model 
described herein for a 2.8 jxm particle suspended in a solution of ionic strength 0.1 

1 0 mM. The legend refers to the ratio of counterion mobility in the double layer to the 
bulk ion mobility. (B) is a plot of the relation between surface potential and 
mobility for the specific relaxation frequency of 1000 Hz. 

Fig. 4 is an illustration of flow fields characteristic of a sandwich geometry 
compatible with LEAPS. 

1 5 Fig. 5 is a flow chart summarizing an analytical procedure for the 

simultaneous determination of particle surface potential and surface conductance 
from experimentally determined values of particle velocity during array assembly 
and particle relaxation frequency. Also indicated are the computer programs 
("Codes") generated to perform the calculations. 

20 Fig. 6 is an illustration of particle movement and of the relative position of 

each particle relative to a position defined as x=0. The result of such particle 
movement is array assembly and each particular particle position represents a stage 
in the assembly process. The particles assemble in thin oxide region of the 
silicon/silicon substrate and the demarcation between thin oxide and thick oxide 

25 defines the position x=0. 

Fig. 7(A) is an illustration of the particle velocity at the interface as a 
function of surface potential, as evaluated numerically on the basis of the theoretical 
model described herein in for several ratios of the counterion mobility to the bulk 
ion mobility, as indicated in the legend. (B) is a plot of the relation between surface 

30 potential and mobility for the specific bead velocity of 0.8 mm/s. 

Fig. 8 is an illustration of the process particle fractionation by selection of 
frequency, co (cdr 1<; go < cor 2 ) from an initially random binary mixture of particles 
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having two different relaxation frequencies (cor 1 and co R 2 , cor 1 > co R 2 ) using a 
substrate surface featuring a half plane of lower impedance adjacent to one of higher 
impedance. 

Fig. 9 is an illustration of the determination of the velocity of particles 
5 experiencing acceleration as they cross a boundary, located at position "0", 
separating regions of high and low electrode impedance. 

Fig, 10 is an illustration of experimentally determined dependence of 
particle relaxation frequency on particle diameter for three different types of 
particles. Experimental details are indicated in the inset to the graphical display. 
10 The results for three different salt concentrations are shown for two ionic diffusion 
coefficients. Results show that the relaxation frequency is not significantly affected 
by a change in salt concentration. 

Fig. 1 1 is an illustration of the process of fractionation of an initially random 
particle mixture composed of larger and smaller particles on the basis of particle 
1 5 size using a substrate surface featuring a strip of lower impedance (bright area) 
surrounded by areas of larger impedance (dark areas). 

Fig. 12 is an illustration of the two successive stages in the process of 
fractionation of an initially random particle mixture composed of particles of 
identical size which display either the protein streptavidin (shown in red) or 
20 displaying a layer of short oligonucleotides attached to the streptavidin coated 
particles by way of a biotin functionality (shown in turquoi se ) . 

Fig. 13 is an illustration of the collection and assembly of HL-60 Human 
Myeloid Cells (-15 |im diameter) in characteristic configurations on a patterned 
silicon chip. 

25 Fig. 14 is an illustration of the collection and assembly of E. coli and staph 

bacterial cells (-1 jum diameter) in characteristic bands on a patterned silicon chip. 

Fig. 15 is an illustration of serial fractionation of a multi-component particle 
mixture (with each particle type possessing a unique relaxation frequenc^ooRVsing 
patterned electrode and sequence of illumination patterns. 
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DETAILED DESCRIPTION OF THE INVENTION 

I. ELECTROSTATIC THEORY OF DIELECTRIC MEDIA 

So that the present invention may be better understood, certain results from 
the electrostatic theory of dielectric media are herein summarized, making general 
reference to Landau, Lifshitz & Pitaevskii, "Electrodynamics of Continuous 
Media", Butterworth-Heinemann, Oxford, UK, 1998. 

LI Ponderomotive Forces 

The free energy of a dielectric medium of dielectric constant 8m filling a gap 
of size H between two parallel capacitor plates of area A carrying a charge Q is 
given by 3o = -ViCofy 2 , where the capacitance, Co = 6m A/H, describes the 
proportionality between applied potential and stored charge, Q= C 0 <|). Under the 
assumption that the electrostatic potential of the plates remains unchanged, the 
introduction of a sphere of radius a « H and dielectric constant e s into the gap 
changes the free energy in accordance with the following expression: 



15 



20 



(1) 



Here, the factor k = (e s - £ M )/(es+2s M ) determines the potential of an 
uncharged dielectric sphere of radius R and dielectric constant 8 S embedded in a 
dielectric medium of dielectric constant 8m and subjected to a homogeneous electric 
field E 0 . 



<f> = -£ 0 rcos 



r>a 



(2) 



In accordance with Eq.(l), as long as 8m< e s , insertion of the sphere into the 
dielectric medium in fact lowers the free energy. Consequently, a sphere placed just 
outside the edges of the plates will be drawn into the space between the plates by a 
force F~ V3. From Eq. (1), with E 0 = <|> IH\ 



25 



F = 2^ M ^V£" 



(3) 
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This ponderomotive force thus arises from the interaction of the polarization 
of the sphere with spatial non-uniformities in the electric field. Eq.(3) shows that the 
force is independent of the sign of the electric field. That is, in a non-uniform 
electric field, a polarizable particle will move either into or out of a region of higher 
5 field, as determined solely by the relative magnitude of the dielectric constants of 
particle and suspending medium. In general, as discussed herein, the polarization 
factor k will depend on frequency: k = £(co). 

1.2 Interfacial Polarization in Electrolyte Solutions 

The polarization of a fluid medium, herein referred to as an electrolyte 

1 0 solution, typically manifests itself in the form of gradients in the spatial distribution 
or orientation of the constituent molecules. Thus, electrolyte solutions are polarized 
near any interface with a bounding (boundary) surface or interface. When such an 
interfacial region is subjected to an external electric field of varying frequency, the 
response of the medium is determined by the processes mediating redistribution or 

1 5 reorientation of the constituent molecules and by characteristic molecular mobilities 
associated with these processes. 

Spatial Distribution of Free Ions: The Double Layer - The response of an 
electrolyte containing mobile charges to the presence of a charged bounding surface 
is the formation of a layer of excess charge adjacent to the surface so as to screen 

20 the field generated by the charges on that surface. Thus, the polarization of the 
interface between a charged particle and an electrolyte containing free ions 
manifests itself in the form of a redistribution of ions so as to form screening shell 
of excess counterfoils. This "double layer" is characterized by a surface potential, <|> s , 
which corresponds to an excess of free charge carriers adjacent to the sphere and by 

25 a characteristic screening length, k" 1 , which is directly related to the ionic strength, 
ViEnjoZi 2 , of the medium; for a 1 :1 electrolyte, k 1 = (s M eol^T/2z 2 e 2 no) ,/2 , where no is 
the molar concentration of the electrolyte and z is the charge on each ion. The local 
excess concentration of ions can affect binding reactions involving molecules 
anchored to the sphere surface. A uniform electric field induces the redistribution of 

30 ions adjacent to and within the "double layer", leading to a phenomenon sometimes 
referred to as "concentration polarization" (Fig 1). The redistribution of the 
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molecules forming the screening layer is determined by translation^ molecular 
mobilities. 

Preferred Orientation of Dipolar Molecules - The response of an electrolyte 
containing dipolar molecules to the presence of a charged bounding surface is the 
5 formation of a layer in which those molecules assume a preferred orientation. As 
with the previous case, the interfacial region of preferred molecular orientation is 
characterized by a surface excess polarization and a characteristic decay length 
which is readily determined from a statistical theory. The local excess orientation 
can affect binding reactions involving molecules anchored to the sphere surface. 
1 0 The redistribution of the molecules forming the screening layer is determined by 
orientational molecular mobilities. In any given electrolyte, excess concentration 
and orientation (as well as higher moments of the molecular charge distribution) 
may contribute to the formation of a screening layer. 

1.3 Frequency Dependence of Field-Induced Polarization 

1 5 An externally applied electric field alters the interfacial polarization of 

electrolytes by inducing the redistribution or reorientation of the constituent 
molecules. Thus, an applied AC electric field induces a time-dependent interfacial 
polarization. As with the recharging of a parallel plate capacitor, the time-dependent 
interfacial polarization corresponds to a transient charging currency ~ d t P. When 

20 the time-dependent interfacial polarization is mediated by the spatial redistribution 
of free ions, the charging current,/ will be composed of a longitudinal and a 
transverse component, respectively reflecting the redistribution of ions within the 
double layer and the exchange of charge between the double layer and the 
surrounding medium. Set forth below are certain results for this case; it is noted that 

25 corresponding results also may be cited for the time-dependent reorientation of 
dipolar molecules. 

The longitudinal component reflects the redistribution of primarily of 
counterions within the double layer in response to the applied field on a 
characteristic time scale i|| - (a + 1/k) 2 /D ~ 1/go r , where a, 1/k and D respectively 

30 denote: the particle radius, the Debye screening length and a characteristic 

diffiisivity of the majority carrier. The characteristic time scale, x = Xjj, is determined 
primarily by the counterion mobility, u - ze/f, where, n and / respectively denote the 
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local counterfoil concentration and drag coefficient. Ionic mobilities in the vicinity 
of the particle surface will reflect the details of the local chemical environment and 
generally will differ significantly from the corresponding value in the bulk medium. 
By invoking the Einstein relation l/f= Dlk^T, the relaxation frequency, (cor, may be 
5 expressed in terms of the counterion mobility, u: 

Dk 2 IF 

■ ; — 1 (4) 

The longitudinal counterion current is balanced by a transverse current 
reflecting the radial transport of co- and counterions across the double layer on a 
time scale x± ~ 1/k 2 D ~ l/a>±. If electroneutrality is to be strictly maintained in the 

1 0 electrolyte adjacent to the double layer, then the exchange of counterions between 
double layer and surrounding electrolyte via radial transport must be accompanied 
by the radial transport of an equal number of coions. The result is a change in the 
local ionic strength in the electrolyte adjacent to the double layer. That is, gradients 
in electrolyte concentration appear on the scale of the particle diameter, and these 

1 5 decay away on a characteristic time scale t ~ 1/cor. As the frequency increases, the 
redistribution of ions within the double layer eventually will no longer suffice to 
attain the fully polarized state; relaxation occurs when co R x - 1 . For a typical 

particle radius of 1 nm, and a typical value co ± - 10 5 /s, the relaxation frequency, 
oo~ co±/(l + Ka) 2 is approximately 1 kHz. 

20 Phenomenological Theory for Polarization Factor k(oo) - Phenomeno logical 
models have been developed to relate the field-induced, frequency-dependent 
polarization, k = k((6) to the particle's surface potential, (j) s , and to the surface 
conductivity, a. The former, for given bulk ionic strength, determines the total 
charge density within the double layer; the latter determines the rate at which the 

25 double layer is able to readjust in response to a time-varying applied field. 

Specifically, a linearized theory, valid for small surface potentials, is available to 
treat the contribution of transverse ionic fluxes to the transient ionic charging 
currents and the charge balance within the double layer to produce an expression for 
the frequency dependence of k = A(co) [Dukhin and Shilov, "Dielectric Phenomena 
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and the Double Layer in Disperse Systems and Polyelectrolytes", Keter Publishing, 
Jerusalem, 1974]. A recent formulation of the linearized theory of the electrokinetic 
transport phenomena underlying concentration polarization, developed for the 
analysis of the low-frequency dielectric response of colloidal suspensions, yields a 
useful analytical form for the frequency dependence of k((o) [Minor et al, J. Colloid 
& Interface Science 206, 397-406 (1998)]. 

Assuming low values of the external field, the potential and the ion 
distribution(s) are linearized around the equilibrium condition existing in the 
absence of an external field. For the region within the suspending medium far from 
the double layer, ion balances (for the i-th species) obey the Nernst-Planck equation: 



Here, the terms preceded by '8' denote the respective first order deviations 
from the corresponding equilibrium values, and f is an effective friction factor 
which governs the ion mobility. Solving Eq. (5) along with the linearized Poisson 
equation, expressions for Srii f ;f and 5<|) f fare found as follows: 



Here, Xa = (1+/) (cox) 1 ; the time scale, x, is set by the characteristic diffusion time 
for transport around the perimeter of the particle: r= a 2 /D =f°a 2 /2k B T 

The polarization parameter, determines the increase in the effective 
dielectric constant of the medium due to the presence of particles and associated 
double layers. To express k = k((6) in terms of system properties of interest, 
boundary conditions are applied to Eq.(6) in the form of a flux balance: tangential 
ion fluxes within the double layer are balanced by radial ion fluxes into and out of 
the double layer so as to ensure that electroneutrality is maintained in the external 
solution at all times. Restricting attention to an electrolyte containing only a single 
coion (superscript "1") and a single counterion species (superscript "2") and 




k h r a 



(5) 




(6) 
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ignoring the coion flux within the double layer as negligible compared to the 
counterion flux, the flux balance at the outer edge ( r = a~) of the double layer is 
expressed in the form: 



(7) 



10 



where V/ denotes a tangential gradient operator. Radial and tangential counterion 
fluxes are related by the constant, Du2, which represents the dimensionless ratio of 
(the counterion-mediated) surface conductivity to (the counterion-mediated) bulk 
conductivity, Du 2 = olac^ a denoting the particle diameter. Du 2 may be expressed 
in terms of the surface potential, <J> S , and a dimensionless counterion mobility, m 2 , as 
follows: 
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AppUcation of the boundary conditions yields the following expression for 
the polarization parameter, &(co): 
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Fig. 2 shows representative plots of £(co) for different values of surface 
potential and for a given ratio of surface ion mobility, w, and bulk ion mobility, u 0 
The point of inflection in each of these plots identifies a peak in the dispersion, 
which in turn identifies the corresponding frequency as the relaxation frequency, 
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cor. The plots in Fig.3 show that for given value of ulu^ © R generally will exhibit a 
dependence on surface potential. That is, a> R may be understood as a function of the 
ionic conductivity, a ~ eun ~ ze 2 nlf, a quantity which in turn depends on both ionic 
mobility as well as local ion density, the latter being controlled by the surface 
5 potential. 

However, Fig, 3 also shows that, for sufficiently large values of the surface 
potential, <j> s , the dependence of cor on surface potential vanishes, yielding instead a 
linear relationship, cor - ulu^ in accordance with the scaling relation, Eq. 4. This 
relation, which formed the basis for the analysis of the data in Fig. 10, thus provides 
1 0 the basis for a direct experimental determination of ionic or molecular mobility from 
a determination of the relaxation frequency. 

H. ELECTROKINETIC PARTICLE TRANSPORT 
AND PARTICLE ARRAY ASSEMBLY 

Light-controlled Electrokinetic Assembly of Particles near Surfaces 
1 5 (LEAPS), as disclosed in W097/40385, enables the formation of planar arrays of 
particles (including on-cue formation) in designated areas of a substrate under real- 
time optical control LEAPS combines the use of electric fields, interfacial 
patterning and/or light to introduce spatial and temporal modulations of the 
impedance of an Electrolyte-Insulator-Semiconductor (EIS) interface. LEAPS 
20 invokes electro hydrodynamic forces which arise in accordance with the lateral 
impedance gradients to control the coupled transport of fluid and suspended 
particles as well as the array assembly process. 

The present invention describes the role of polarization-induced forces that 
contribute to particle transport and array assembly under conditions characteristic of 
25 LEAPS. These models guide the selection of conditions for the fractionation of 
particle mixtures in response to combinations of time-varying electric fields. These 
models also provide methods for the analysis of coupled fluid -particle transport in 
order to determine characteristic particle properties of interest. 

in. FLUID FLOW, PARTICLE COLLECTION AND 
30 ARRAY ASSEMBLY 

Light-controlled Electrokinetic Assembly of Particles near Surfaces 
(LEAPS) is described in PCT International Application No. PCT/US97/08159 
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(published as WO 97/40385) and U.S. Serial No. 09/171,550, filed October 26, 
1998 (issued as U.S. Patent 6,251,691 on June 26, 2001) which are incorporated by 
reference in their entirety. Also incorporated herein by reference in their entirety are 
U.S. Serial No. 09/397,793, filed September 17, 1999 (published in WO 01/20593) 
5 (describing system and method for programmable illumination pattern generation), 
as well as U.S. Serial No. 09/320,274, filed May 28, 1999, and PCT International 
Application No. PCT/US00/14957 (published as WO 00/73777) (relating to array 
cytometry). LEAPS involves the interplay between several electrokinetic forces to 
induce fluid flow, particle transport and array assembly. For example, in a sandwich 

10 geometry, an AC electric field applied across a fluid gap generates rectified 

electroosmotic flow in the vicinity of lateral gradients of the Electro lyte-Insulator- 
Semiconductor (EIS) impedance that are set up by interfacial patterning or spatially 
varying illumination (see Fig. 4). Generally, these flows are observed under the 
condition <o < toj_, the latter characteristic frequency being determined by molecular 

1 5 processes governing the response of the polarization of the EIS interface. Under the 
action of these flow fields, particles are "collected" over macroscopic distances and 
are swept from regions of higher toward regions of lower interfacial impedance by 
viscous drag. 

Polarizable particles, while being swept along by the electrokinetic flows, 
20 simultaneously respond to the applied electric field; specifically, if particles are 
charged, the applied field induces local charging currents which in turn create a 
substantial concentration polarization in accordance with the mechanism previously 
discussed herein. This field-induced polarization determines the effective dielectric 
constant of polarizable particles and so dominates the ponderomotive forces exerted 
25 on such particles in regions of rapid spatial variations in electric fields such as those 
produced near steps or other steep gradients in interfacial impedance. 

Accordingly, particles approaching a step in EIS impedance come under the 
influence of the gradient V(E-E) and accelerate across the step into the low 
impedance region where they are slowed, e.g., by the drag force exerted by the 
30 fluid. The particles enter into regions of low impedance and form arrays when 

0) < cor ("aggregation"); in contrast, when co > cor, particles accumulate on the high 
impedance side in the vicinity of the boundary. Illumination also can induce 
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gradients in impedance and, provided that co < co R , array formation in illuminated 
areas of the electrode. 

In certain embodiments of the present invention, the frequency dependence 
of the array formation process permits the experimental determination of the 
characteristic frequency, co R . For a homogeneous population of particles, the 
determination of the relaxation frequency is determined by inducing a transition in 
the array configuration: particle arrays assembled under the condition co < co R 
assume a close-packed configuration; when the frequency is then timed to co > co R , 
arrays undergo a transition to an expanded configuration [Chau et al, Proc. SPIE, 
3877, 36-42] . Accordingly, the frequency of the applied AC electric field is first set 
to a sufficiently low value so as to enable the assembly of an array within a 
designated area and is then set to progressively higher values to trigger the transition 
to an expanded array configuration. 

For a heterogeneous mixture of particles disposed on a patterned electrode 
surface, the entire set of relaxation frequencies characterizing the types of given 
particles maybe determined as follows. The frequency of the applied AC electric 
field is first set to a sufficiently high value so as to prevent transport of any particles 
into designated areas of low impedance, and is then set to progressively lower 
values to permit array assembly for those particles satisfying the criterion co < cor. 

Particle Size - Using this transition in array configuration to determine co R , 
the relation between co R and particle size may be determined for different electrolyte 
concentrations and for particles displaying different surface properties (see Example 
2). 

Surface-Chemical Composition - For given size, cor also varies with particle 
surface chemistry. For example, when proteins are coupled to bead surfaces, co R 
decreases; when oligonucleotides are coupled to protein-coated beads, or, increases. 
In general, both the local density of charge carriers, n, as well as their ionic 
mobility, w, maybe affected. 

Non-aqueous Electrolytes - Rectified field-induced fluid flow as well as 
particle transport and array assembly, the latter accompanied by an acceleration 
across impedance steps, also occur in non-aqueous suspensions. The relaxation 
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frequency of particles assembled in DMSO and acetonitrile is shifted to lower 
values when compared to aqueous suspensions. 

IV. FRACTIONATION OF HETEROGENEOUS PARTICLE 
MIXTURES 

5 The characteristic dependence of co R on particle size and surface chemical 

composition, manifesting itself in its effect on particle transport and array assembly 
under conditions characteristic of LEAPS as described herein, provides the basis for 
an on substrate (e.g., on-chip) method to fractionate complex heterogeneous 
mixtures of particles including suspensions of cells. The underlying principle of 

10 operation of the method of the present invention is to form designated areas of low 
impedance of an EIS interface, either by interfacial patterning or by illumination as 
provided in LEAPS, and to select a frequency, or sequentially select multiple 
frequencies, of the applied voltage, under conditions compatible with the 
aggregation of at least one constituent subset among a given, heterogeneous set of 

15 particles. 

To fractionate a mixture composed of two different types of beads that differ 
in their respective relaxation frequencies, the frequency of the applied voltage is set 
to an intermediate value to selectively assemble arrays of beads of the higher 
relaxation frequency. In accordance with the dependence of relaxation frequency on 

20 particle size, as disclosed herein, the condition co R L < <o < co R s ensures that only the 
smaller particles (characterized by a frequency cor S ) enter the region of lower 
impedance, while larger particles (characterized by a frequency cdr L ) align on the 
high impedance side of the feature (see Example 4). Similarly, the dependence of 
the relaxation frequency on surface chemical composition can be used to fractionate 

25 particles (see Example 5). 

The binding of ligands to the surface of functionalized beads that produces a 
change in surface chemical composition can be detected on the basis of the 
corresponding shift in relaxation frequency. In its simplest format, the process is 
illustrated by the selective aggregation in a low impedance area of the electrode of 

30 only one of an originally randomly mixed set of two types of particles, with 

respective relaxation frequencies gor 1 < a> R 2 , by selection of the operating frequency, 
(Or < ©r < (Or 2 , as illustrated in Fig. 8. 
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This fractionation process of the present invention may be generalized to 
multi-component mixtures by progressively selecting frequencies to isolate and 
assemble one constituent subset of beads. The method of the present invention also 
applies to cells whose differential adhesion properties provide an additional 
5 discriminating feature for LEAPS-mediated fractionatioa By the method of the 
present invention, different types of bead and cells are fractionated and sorted 
without the need for complex and expensive equipment such as flow cytometers. 

As with the method disclosed herein, polarization-dependent forces also play 
a role in dielectrophoretic methods for fractionation of heterogeneous cell 
1 0 populations. However, the methods of present invention provide several distinct 
advantages, including: 

- the combination of electrokinetic flow and polarization- 
dependent particle responses in a collective process of particle array assembly, 
thereby providing for the integration of particle transport with fractionation in 

1 5 designated areas of an electrode; 

- fractionation on the basis of surface chemical properties, thereby 
providing for the detection of changes in the surface composition of particles or 
cells in the course of binding assays; 

- sequential multi-step fractionation within substantially the same 
20 area of a substrate under optical monitoring; 

- fractionation by multiple parameters including size, surface 
composition and selective adhesion; 

- multi-dimensional fractionation, e.g., 

frequency dispersion of array assembly, followed by selective adhesion; 
25 - parallel, multi-contrast mode of detection via optical imaging; 

- ease of use due to simple handling of electric contacts and; 
flexible electrode patterning by standard methods of semiconductor processing; 

- compatibility with dielectrophoretic separation processes; 

- optical programmability as provided by Seul and Chau [Chau et 
30 al, Proc. SPIE, 3877, 36-42]. 

Several related distinctions may be drawn between the method of the present 
invention and dielectrophoretic methods of particle separation including: 
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Underlying Process - The concentration polarization process invoked by the 
present invention, detectable by dielectric spectroscopy of particle suspensions 
[Minor et al, J. Colloid & Interface Science 206,397-406 (1998)], probes molecular 
configurations and transport phenomena associated with particle surfaces. For 
5 particles within the range of sizes of typical interest to the applications considered 
herein (0.1 |nm - 10 jam), this process is associated with characteristic relaxation 
frequencies in the range of a few hundred Hz to a few tens of kHz, a range 
comprising AC field-driven electrokinetic flow phenomena associated with the 
electrode surface [Yeh et al, Nature 386, 57-59 (1997)]. 

10 In contrast, representative implementations of dielectrophoretic methods for 

fractionation of particle mixtures developed by Pethig's group [G. H. Markx, M. S. 
Talary, R. Pethig, J. Biotechn 32, 29-37 (1994)] and by Gascoyne's group [Y. 
Huang, X-B. Wang, F. F. Becker & P. R. C. Gascoyne, Biophys. J. 73, 1 1 1 8-1 129 
(1997)] use the frequency-dependent transition from a positive to a negative 

1 5 dielectrophoretic response of cells at high frequency (typically exceeding 100kHz) to 
distinguish between cell types; this transition, detectable by the response of 
individual cells to rotating electric fields in electrorotation measurements, reflects 
the relative contributions of the dielectric permittivity and membrane conductance 
to the overall dielectric response of cells [Jones, "The Electromechanics of 

20 particles", Cambridge University Press, Cambridge, UK 1995]. In contrast to the 

collection of particles over macroscopic distances afforded by electric-field-induced 
fluid flow in accordance with the present invention, the prior art relating to 
dielectrophoretic methods of particle fractionation does not invoke electric field- 
induced fluid flow in the fractionation process. 

25 Electric Field Configuration - As with LEAPS, the configuration of electric 

fields according to the method of the present invention is produced by applying an 
AC voltage normal to the EIS interface, for example by adopting a sandwich 
configuration in which the voltage is applied between a pair of parallel planar 
electrodes confining the electrolyte solution. That is, the applied field is nominally 

30 normal to the electrode surfaces. In contrast to the methods of the present invention, 
voltages in these dielectrophoretic methods of separation are not applied across the 
fluid gap, but are instead applied between pairs of electrodes deposited on the same 
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surface. Because the forces of interest in these methods of the prior art are of very 
short range, arising primarily from the spatial inhomogeneities of the "fringe" fields 
at the electrode edges, designs typically employ sets of closely spaced metal 
electrodes of opposite polarities. 
5 Pethig and colleagues use interdigitated, "castellated" microelectrode 

designs and scan frequency to separate cells by positive and negative 
dielectrophoresis; for example, a range of frequency can often be found in which 
viable cells experience positive dielectrophoresis while non-viable cells experience 
negative dielectrophoresis. Differential forces of retention strongly confine cells 

1 0 under positive dielectrophoresis, permitting cells experiencing negative 

dielectrophoresis to be swept away by an externally provided lateral fluid flow. 
Gascoyne and colleagues have devised a dielectrophoretic levitation metho v d in 
which traveling wave dielectrophoresis is combined with externally driven laminar 
flow along the electrode surface in a field-flow fractionation process in which the 

1 5 differential dielectrophoretic lift of different types of cells exposes them to different 
flow velocities of a parabolic flow profile maintained in the separation chamber. 
Entrained by the flow, cells are detected, one at a time, as they traverse a 
downstream cell counter. 

Simultaneous Transport and Assembly - The LEAPS-enabled fractionation 

20 method of the present invention combines electrokinetic flows with concentration 
polarization in a collective process to integrate flow-mediated particle transport with 
the differential response to ponderomotive forces to bring about both the long-range 
transport of particles ("collection") as well as fractionation in the process of array 
assembly. The frequency ranges typically employed in dielectrophoretic analysis 

25 exceed the characteristic relaxation frequencies associated with electric field- 
induced electrokinetic flows, a fact which has profound consequences for the design 
of electrode configurations permitting the short-ranged dielectrophoretic forces to 
be utilized for particle collection. Voltage switching strategies (as for example in 
traveling wave dielectrophoresis), enabling serial transport and fractionation of 

30 particles, however, suffer from the disadvantage that it requires complicated 
electrode configurations and design. 
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Optical Programmability - The method of the present invention is unique in 
providing for real-time control of fractionation processes by means of light-induced 
interfacial impedance modulations afforded by the programmable generation of 
illumination patterns. 

V. PARTICLE ANALYSIS 

The present invention includes a method of analyzing particles based on the 
analysis of particle transport in the course of particle array assembly in accordance 
with LEAPS. To that end, the equation of motion of a spherical particle including 
viscous drag and ponderomotive forces may be written in the form: 



where Vb and Vf respectively denote the velocity of the particle and the fluid medium 
and K((6) = 27ia 3 e M &, where k is defined in connection with Eqs (1) and (2). The 
ponderomotive force will vanish as k -» 0 or 8s -> Sm. Reference to Fig. 2 shows 
that this condition corresponds to relaxation: as the frequency of the external electric 
field approaches co R , the fully polarizable state is no longer attainable. In general, 
v max will thus depend on surface potential, ion mobility and frequency. 

This equation may be used to analyze the motion of particles in the course of 
particle collection and array assembly by implementing the following steps: 

1 . For a given EIS patterning geometry, specified, for example, 
in the form of the electrode surface potential, <|)s(x, y) = <J)(x, y, z=0), solve the 
Laplace equation to obtain the potential and the electric field distribution, E (x, y, z) 
= -V<j)(x, y, z); for several relevant configurations, analytical models of <t> s (x, y) = 
(|)(x, y, z=0) are available (see Example 1); 



the slipping plane of the electrode; this serves as the boundary condition for the 
evaluation of the flow field above the electrode, further assuming the fluid to be 
static at infinite distance from the electrode; 




3. 



Evaluate the term V(E-E); 

Evaluate the electrokinetic fluid velocity, of Vf 1 ** = -< 



e/uV((), at 
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4. Determine K(g)) using the expression: K(<i)) =2to 3 8 M (s s (o))- 
eivi)/(es(co)+2£M); 

5. Determine the particle velocity, dxddt = Vb(t) + Vf(t), in the 
particle moving frame; the maximum (relative) particle velocity, Vmax> is attained in 
the vicinity of steep impedance gradients; determine the particle position, x*> = Xb(t); 

6. Determine (using the functions obtained in steps 4 and 5), the 
particle velocity as a function of position: dx\Jdt =J(xb). 

In certain embodiments of the present invention, the algorithm detailed in 
the flow chart presented in Fig. 5 is used to convert the experimentally determined 
values for the relaxation frequency, <or, and the frequency-dependent maximum 
particle velocity, Vmax, attained near impedance gradients (see e.g. Fig. 6) into the 
particles' surface ("zeta") potential, <(> s , and surface conductivity, a. As indicated in 
Figs. 3 and 7, both cor and v max depend on surface potential and ionic mobility and 
the determination of a pair of parameters (<|> s , u) such that co R = co R exp and 
Vmax - Vmax exp is carried out using numerical analysis as summarized in the flow 
chart. The present invention also includes a method for particle analysis that allows 
simultaneously determination of morphological parameters such as size and shape 
along with the electrical parameters in a miniaturized imaging format. 

This method provides the basis for a quality control technology that is 
particularly sensitive to surface fixnctionalization. In addition, the method provides 
the basis to detect the binding of biomolecules to "matching" molecules displayed 
on the bead surface. That is, the change in surface chemical composition that results 
from the binding of certain molecules to particle or cell surfaces is detected without 
the need for fluorescent or other optical labeling. 

The present invention will be better understood from the Experimental 
Details that follow. However, one skilled in the art will readily appreciate that the 
specific methods and results discussed are merely illustrative of the invention as 
described in the claims which follow thereafter. 
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EXAMPLES 

General Experimental Conditions: 

Preferred particle array assembly is carried out in fluidic microcells which 
are formed by sandwiching a double-sided Kapton spacer of -100 jxm thickness 
5 between a 1 cm x 1 cm substrate that can be patterned electrically, such as a silicon 
chip (n-type, capped either by a uniform or a lithographically patterned thin SiC>2 
layer), also serving as the bottom electrode, and a glass cover slip coated with 
indium tin oxide (ITO) to a typical sheet resistance of 1400 Ohm Square serving as 
the top electrode. AC voltages of 1 -10 Vp. p , in a frequency range from 100's of hertz 

10 to a few kilohertz are applied between the electrodes across the fluid gap. Particle 
transport and assembly were monitored by video microscopy permitting frame 
capture and digitization of frames for further analysis. 

Aqueous Suspensions - Prior to assembly, particles stored in buffer are 
washed and resuspended in deionized and ultrafiltered water (conductivity < 50 S 

15 cm" 1 ); in some cases, salt (< ImM) is added to provide buffering. Anionic and 

cationic particles from 0.5 jxm to 10 |im in diameter, composed of silica or modified 
polystyrene and functionaUzed with a variety of chemical surface groups, as well as 
functionalized core-shell particles obtained from a variety of manufacturers have 
been used. These particles were further modified to display oligonucleotides, DNA 

20 fragments or proteins. 

Non-Aqueous Suspensions - Illustrative of these systems are: 
DMSO(8 DM so=47.24): 4.4 jam carboxylated PS particles (Dyno, Oslo, Norway); 2.1 
^m carboxylated PS particles stained with fluorescein (Molecular Probes, Eugene, 
OR); in both these cases, stock particles originally suspended in water/buffer, were 

25 resuspended in DMSO by way of several successive cycles of centrifugation and 
resuspension to exclude water. Glass beads (3-10 jim in diameter) (Polyscience, 
PA); stock particles, supplied in dry powder form, were resuspended in the solvent 
of choice prior to assembly. Acetonitrile (s a cetomtriie 36.64): 1 jim nonporous oxirane 
acrylate beads (Sigma, MO); procedure analogous as in the case of glass beads. 
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Example 1- Particle Array Assembly: Half-space Geometry 

Fig. 6 displays bead velocities measured in the course of array assembly of 
oligo-dT(25) particles of diameter 2.8 ^m (obtained from Dynal, Oslo, Norway) 
using kg s 45 (see also the conditions shown in the inset of the lower panel) and 
5 invoking video microscopy. Analysis of the "streaking" in the snapshots in the 
upper portion of Fig. 6, taken with 1/30 s exposure, reveals the significant 
acceleration of particles traversing the impedance step, created by a variation in 
oxide thickness from -1000 A to -100 A. Visual inspection also reveals that 
particles are simultaneously displaced toward the substrate. The instantaneous 

10 velocities extracted by image analysis of a series of snapshots are plotted in the 
lower panel of Fig. 6. Arrays of assembled particles act in the manner of low 
impedance regions of the substrate: they generate in-flow directed toward the array 
perimeter and attract additional particles which accelerate toward the array. 

The numerical solution of the equation of motion, given an impedance step 

15 in a "half-space" geometry, reveals a strong dependence of v max on surface potential 
for several ratios of u/u«> and for given particle size, as illustrated in Fig. 7. This is 
as expected: a high surface potential corresponds to a large induced dipole moment 
and hence to a large ponderomotive force; in contrast, the dependence on ion 
mobility is weak. 

20 The equation of motion takes into account the ponderomotive forces, the 

drag forces and also the moving frame of reference of the fluid. The path of any 
bead can be simulated given the experimental conditions such as voltage, gap size, 
ionic strength, ion mobilities, oxide thicknesses and given pertinent bead properties 
such as size and surface potential. In general, it is desirable to apply the numerical 

25 model to fit experimental data and so to improve the accuracy of the determination 
of v max . 

A typical solution of the equation of motion, generated numerically by 
representing the "half-space" geometry by a surface potential profile in the form 
<|> s = A tanh(icc), is shown in Fig. 9. The origin of the graph represents the location of 
30 the impedance step. The position of the bead is negative when it is in the thick oxide 
region, where the velocity of the bead is determined by the fluid velocity at any 
given point. However, when approaching the impedance step to within a distance of 
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the order of the Debye length, the bead experiences the aforementioned 
ponderomotive force which causes it to accelerate across the step relative to the 
fluid. Given the short range of ponderomotive forces (~Vi? 2 ), the acceleration ceases 
within a short distance of the step, and the bead experiences a drag force which 
5 decelerates it to the velocity of the surrounding fluid. At larger distances from the 
step, the lateral fluid motion ceases and consequently, the bead comes to rest unless 
the ponderomotive forces surrounding pre-existing bead arrays or clusters provide 
an additional acceleration toward the array boundary. 

Example 2 -Dependence of Relaxation Frequency, cor, on Particle Size and 
10 Composition 

Fig. 10 displays the results of determining the relaxation frequency by 

inducing an array expansion for a variety of particle sizes and for different types of 

particle surface chemistries. In a first group, particles included were carboxylate- 

modified or amine-modified polystryrene particles (obtained from Dyno, Oslo, NO) 

1 5 as well as "oligo-dT(25)" particles (obtained from Dynal, Oslo, NO). In a second 
group, particles were coated with a monoclonal antibody against interleukin-6 
(obtained from R&D Systems, Minneapolis, MN). The anticipated scaling behavior 
of the relaxation frequency, o>r 1/2 ~a (particle size), as discussed herein, is 
confirmed by the experiments. 

20 The solid lines represent a fit of the data to a simple expression discussed 

herein and produces a diffusivity parameter, D, which reflects the local ion- 
mediated surface conductance, as discussed herein. Invoking Eq. 4, the data in Fig. 
10 may be interpreted to indicate that, for given a and k, the binding of proteins, or 
other large molecules to the bead surface reduces the surface conductivity, a ~ e 2 n/f, 

25 introduced herein. 

Example 3 - Particle Analysis 

From experiments illustrated in Figs. 6 and 10, values for the maximum 
velocity of beads traversing the step from thick oxide to thin oxide, 
v max mm/s, and the relaxation frequency, a> R - 2 kHz, were determined. The 
30 theoretical model and algorithm underlying Figs. 3, 5 and 7, yield corresponding 

values of <|> s - 0.5V and u/u^ - 3.9 x 10°; this combination in turn yields a value of 
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a a - 0.15 S/m and Du - 290; here, u„ -0.02 S m 2 /C, as determined from the 
tabulated molar conductivity for Na + . The value for § s may be compared to the 
surface potential of -330 mV corresponding to the titratable charge of carboxylate- 
modified polystyrene particles of 1 COO 74 A 2 or 400 jaC/cm 2 ; the results obtained 
5 indicate that the effective surface charge of oligonucleo tide-modified particles 
substantially exceeds that of carboxylate-modified particles. The value for cor 
indicates the surface mobility of small cations relevant to the experiments to be 
smaller by a factor of > 100 than the mobility in the bulk; however, the surface 
conductivity still exceeds the bulk conductivity by two orders of magnitude, 
1 0 reflecting the high value of the surface potential. 

Example 4 - Fractionation of Particle Mixture Using Size-Dependence of (Or 

Fig. 1 1 illustrates the fractionation of a binary mixture of smaller (2.1 \im 
diameter, characteristic frequency or S ) and larger (4.5 jum diameter, characteristic 
frequency a> R L ) particles using a surface design in which a region of lower 

1 5 interfacial impedance is created by a variation in oxide thickness from -1000 A to 
-100 A within the stripe indicated by lighter shading (see also Example 1). When 
the frequency, co, of the applied voltage adjusted such that G) R L < co< co r s , only the 
smaller particles experience an acceleration toward the low impedance region where 
they assemble into an array. By narrowing or widening the stripe, flow along the 

20 stripe also may be generated and may serve as the basis for two-dimensional 
fractionation in the cross-flow so realized. 

Example 5 - Fractionation of Particle Mixture Using Composition-Dependence 
of (Or 

Fig. 1 2 illustrates the fractionation of a binary mixture of red and blue 
25 particles of identical size but differing surface chemical composition. The relaxation 
frequency, a)R bare , of carboxylate-modified particles, tagged with a green dye, is 
lowered by binding of protein, tagged with a red dye, to a characteristic frequency, 
©R mod . Under the condition co R mod < (o < ©R bare , only the unmodified particles enter 
the stripe region of lower impedance, while modified, protein-coated particles align 
30 on the high impedance sides of the feature. The magnitude of the shift in 
characteristic frequency is determined by the amount of bound protein. 

689616 vl 



-30- 

Example 6 - Detection of Protein Binding to Beads or Cells without the Use of 
Tags 

A molecular binding event on a known particle surface causing a change in 
the surface chemical property (i.e. a change in the dielectric property of the particle 
5 plus bound molecule over that of the particle) can be easily be detected utilizing the 
fact that the effect of the binding is reflected in change in the relaxation frequency 
of the particle. For example, Figure 1 0 illustrates the experimentally evaluated 
change in relaxation frequency of a particle which display, on its surface, either 
carboxylate groups or carboxylate plus captured layer of antibody. The binding of 

10 antibody to the particle surface lowered the characteristic frequency, oc>mod 9 of 
modified particles relative to the characteristic frequency, ©bare, of unmodified 
particles. Similarly, the binding of biotinylated oligonucleotide probes on the 
surface of a streptavidin functionalized particle may be detected. The binding of the 
oligonucleotide fragment increased the characteristic frequency, aw, of oligo 

1 5 bound particles relative to the characteristic frequency, ow, of unbound particles. 

Example 7 - Frequency-Dependent Collection and Assembly of Human 
Myeloid Cells 

Fig. 13 illustrates the collection and enrichment of HL-60 human myeloid 
cells (CCL-240, obtained from ATCC and maintained for experiments in 8.5% 

20 sucrose/0.3% dextrose) in designated areas of a patterned silicon substrate. Using 
silicon chips with patterned oxide surfaces featuring a stripe of thin oxide (< 100 A) 
typically 100 -200 ^m in width delimited by "steps" in impedance indicated by the 
change in brightness, application of a voltage of 5V peak-to-peak, with frequency 
set to the values indicated in the two panels in the figure, produces aggregates of 

25 either adjacent to, or within the designated stripe, reflecting the selection of the 
frequency in the applied voltage. 

Example 8 - Assembly of Bacterial Cells and Fractionation by Selective 
Adhesion 

Fig. 14 illustrates the collection and enrichment of E.coli and 
30 Staphylococcus epidermis bacterial cells, serving as respective model systems for 
Gram-negative and Gram-positive bacteria, in different regions of a patterned 
substrate (here a silicon chip). That is, using chips with patterned oxide surfaces 
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featuring a strip of thin oxide (< 100 A) typically 100 -200 \im in width, 
application of a voltage of 5 V peak-to-peak across a fluid gap of 100 jam and 
frequency set to 2 kHz produces "bands" of cells in positions reflecting the 
preferential adhesion of S. epidermis near the "step" in impedance indicated by the 
5 change in brightness and E.coli in the center of the strip. The position of the bands is 
thus characteristic of bacterial type. 

Example 9 - Strength of Adhesion 

Cell rolling along a surface is the initial step in many kinds of inflammatory 
response. The present invention provides for the experimental determination of the 

1 0 strength of adhesion between cells of interest or ligand-functionalized 

microparticles disposed on a receptor-functionalized electrode surface. For example, 
given a patterned electrode surface as described, low-impedance regions can coated 
with E-Selectin to investigate the assembly of HL-60 cells (see example 7). Under 
the action of LEAPS-mediated flow, cells, while being transported toward low- 

1 5 impedance regions of the electrode surface, exhibit rolling and sticking, the latter 
representing a manifestation of a positive adhesive interaction between cells and 
surface. In-vitro assays have been described which use the rolling of cells or ligand- 
functionalized microparticles on a receptor-functionalized surface that is exposed to 
tangential fluid flow to estimate the strength of adhesion between cells or 

20 microparticles and the surface [Brunk, D.K., and Hammer, D.A., Biophysical 

Journal, 72,1997, 2820-2833; Hammer, D.A., and Apte, S.M, Biophysical Journal, 
63,1992, 35-57; Springer, T. , Cell, 76, 301, 1994.]. The analysis of rolling and 
sticking of HL-60 cells on E-Selectin coated glass wafers under the presence of 
pressure driven flow has been used for anti-inflammatory drug screening [Mossavi 

25 et al., Proceedings of MicroTAS ' 98 Workshop, D. J. Harrison and A. van den Berg 
(Eds.), 69-72, 1998]. The present invention provides a superior mode of flow 
control by way of LEAPS. For example, a receptor-functionalized electrode on 
which cells or microparticles of interest have been disposed in predesignated areas 
can be exposed to LEAPS-mediated fluid flow (see Fig. 4) whose direction and 

30 strength are readily controlled. This permits the evaluation of the correlation 

between the strength of the flow-induced shear stress acting on the adsorbed cells or 
particles and the efficiency of the hydrodynamic detachment, which in turn permits 
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an estimation of the adhesive forces operative between the cell or particle and said 
surface. 

Example 10 - Multi-dimensional Fractionation 

In multi-dimensional fractionation, a mixture of particles is separated into 
5 subpopulations in two or more stages in accordance with multiple discriminating 
parameters. Here, in addition to the frequency dispersion of field-induced 
polarization other attributes such as differential selective adhesion etc. can be used. 
For example, a mixture of HL-60, E. coli and S. epidermis can be fractionated using 
patterned chips featuring a strip of thin oxide (< 100 A) typically 100 -200 fim in 
1 0 width. Application of a voltage of 5 V peak-to-peak across a fluid gap of 1 00 |im 
and frequency set to 2 kHz produces "bands" of cells in positions reflecting the 
preferential adhesion of S. epidermis near the "step" in impedance and E.coli in the 
center of the strip. Under this condition the HL-60 cells form linear chains at the 
thick-thin oxide boundary which extend to the thick oxide region ( see example 7). 

1 5 Example 11- Sequential Light-directed Fractionation of Multi-component 
Mixtures 

A multi-component particle mixture in which each particle population 
displays a unique relaxation frequency (cor 1 ) may be fractionated using a 
combination of a patterned electrode and illumination as illustrated in Figure 1 5. 

20 Given relaxation frequencies (cor 1 ) in the range coi ow < ©r 1 < cohigh, the frequency of 
the applied AC electric field is scanned from the high limit (cohigh) to the low limit 
(coiow) of the range; particles are sequentially swept into the thin oxide/low 
impedance collection region starting with the ones with higher relaxation 
frequencies. Once inside, each population is guided by light into separate thin oxide 

25 pads which are connected to the initial collection zone via narrow channels. This 
particular patterning will confine beads in the respective pads. 

Example 12 - Encoding of Particle Arrays using Engineered Particle 
Polarizability 

Dielectric "tagging" to produce particles of desired polarizability and 
30 corresponding relaxation frequency may be used to produce a specific dependence 
of the particle response to LEAPS. Dielectric tagging may be accomplished by way 
of introducing labels to the exterior or interior of the particle. For example, 
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metallization of particles changes their polarizability. This approach may be 
extended to the engineering of a particle whose polarizability is light dependent. For 
example, using light-sensitive tags such as TiCh nanoparticles, the polarization of a 
Ti0 2 -tagged particle is increased exclusively when illuminated with light [Y. 
5 Komoda, Tata N. Rao, and A. Fujishima, Photoelectrorheology of Ti0 2 nanoparticle 
suspensions. Langmuir, 1997, 13, 1371-1373]. 

Although the invention has been described above with reference to examples 
and to preferred embodiments, it will be appreciated that the invention may be 
embodied in other forms or carried out in other ways without departing from the 
10 spirit or essential characteristics thereof. The above description is therefore to be 
considered in all respects illustrative and not restrictive, and all changes which come 
within the meaning and range of equivalency are intended to be embraced therein. 
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